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Abstract 



Metabolic control analysis (Kacser & Burns (1973). Symp. Soc. Exp. Biol. 27, 65-104; 
Heinrich & Rapoport (1974). Eur. J. Biochem. 42, 89-95) was developed for the under- 
standing of multi-enzyme networks. At the core of this approach is the flux summation 
theorem. This theorem implies that there is an invariant relationship between the control 
coefficients of enzymes in a pathway. One of the main conclusions that has been derived 
from the summation theorem is that phenotypic robustness to mutation ( e.g. dominance 
) is an inherent property of metabolic systems and hence does not require an evolution- 
ary explanation (Kacser & Burns (1981). Genetics. 97, 639-666; Porteous (1996). /. 
theor. Biol. 182, 223-232 ). Here we show that for mutations involving discrete changes 
( of any magnitude ) in enzyme concentration the flux summation theorem does not hold. 
The scenarios we examine are two-enzyme pathways with a diffusion barrier, two en- 
zyme pathways that allow for enzyme saturation and two enzyme pathways that have 
both saturable enzymes and a diffusion barrier. Our results are extendable to sequential 
pathways with any number of enzymes. The fact that the flux summation theorem cannot 
hold in sequential pathways casts serious doubts on the claim that robustness with respect 
to mutations is an inherent property of metabolic systems. 
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1. Introduction 



How do changes in enzyme properties alter the physiological phenotype? Such knowl- 
edge is a key component in understanding the relation between genotype and phenotype. 
The two main approaches to this problem have developed into fields known as metabolic 
control an alysis ( MCA ) (|Kacser & Burns, 1973|; Heinrich & Rapoport, 1974 ; Kacser 
et ai, 1995 ) and biochemical systems theory ( BST ) (Savageau, 1976 ; Savageau & Sor- 



ribas, 1989). Results from MCA have had extensive influence in biochemistry, genetics 
and evolution ( for an overview see the entire issue of J. theor. Biol, 182:3 , 1996 ). The 
cornerstone of the MCA approach is a theoretical result referred to as the flux summa- 
tion theorem. The biological interpretation of the summation theorem is that there are 
systemic constraints inherent in metabolic pathways; these constraints limit the magni- 
tude of the effects that changes in enzyme activity can have on flux through a pathway. 
The existence of such constraints has two implications. In the first place, the summa- 
tion theorem implies that in general the control of flux in a pathway is shared between 



enzymes. Hence rate limiting enzymes are rare flKacser, 1995[ ). Secondly, the theorem 
implies that on the average, mutations that change enzyme concentrations will have a 
small effect on flux. This implication has important consequences in evolutionary theory 
and genetics. It means that phenomena that fall under the rubric of phenotypic robust- 
ness to mutations ( such as selective neutrality, canalization and dominance ) are inherent 
properties of metabolic systems and are not a result of evolution. This MCA argument 
was originally made in reference to the case of dominance, stating that dominance is an 
inhe rent property of metabolic pathways and not a result of evolution (K acser & Burns, 



1981; |Keighfley, 1996| ; [Porteous, 1996[ ). 

At the core of the MCA approach is a measure of the control exerted by each enzyme 
on flux through a pathway. The control coefficient Cf measures how important each 
enzyme is in its ability to affect steady-state flux via changes in enzyme concentration. 



In its original formulation (Kacser & Burns, 1973) the control coefficient was defined as 



j_5J 6Ei 



(i) 



where J is the steady state flux of metabolites through the pathway (net rate of product 
formation at the end of a pathway), E\ is the concentration of enzyme i, 5Ei is a discrete 
change in enzyme concentration and 5 J is the resultant change in flux. Hence Cf is 
basically a ratio of percentage change in a pathway's steady state flux to the corresponding 
percentage change in enzyme concentration. 

Theoretical results that contradict or object to different conclusions derived within 



MCA have been posed several times (|Cornish-Bowden, 1987|; piersch, 1988[ Savageau 
& Sorribas, 1989; [Savageau, 1992fc Prossniklaus et al, 1 996] ; [Kholodenko etal, 1998[ ). 
In each case the possible misrepresentation of non-linearities in MCA has appeared in 
a different guise. In the first place, Cornish-Bowden (1987) showed that in a sequential 
pathway if the maximal rate V max of consecutive enzymes are sequentially decreasing 
then dominance is not a necessary property of the pathway. Hence the possibility exists 
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that dominance can evolve. This objection was rejected in MCA on grounds that the 
specific arrangement of kinetic values suggested is a special case that is very unlikely 
to occur by chance ( Kacser, 1987 ). A second set of objections was then put forth by 
Savageau and Sorribas (1989, 1992). They argued that in pathways exhibiting non-linear 
behavior that arises from properties such as enzyme-enzyme interactions, feedback loops 
or non-sequential pathway structure it can be shown that dominance is not inherent to the 
pathway. Based on disagreements on the mathematical models used, this second objec- 
tion was rejected by MCA proponents ( Kacser, 1991 ; Kacser, 1995 ). Nevertheless more 
recent wor k indirectly corroborates this objection (|Grossniklaus et ai, 1996|; K holodenko 
etal, 1998). 

The present paper deals with problems involving the biological in-applicability of 
MCA precepts. We show that for discrete changes ( of any magnitude ) in enzyme con- 
centration, the flux summation theorem does not hold. This means that contrary to the 
MCA allegations, there are no a-priori constraints that would require the magnitude of 
mutational effects to be small. Furthermore these results re-open the possibility that 
phenotypic robustness to mutations ( e.g. dominance ) is not an inherent property of 
metabolic pathways and can be a result of evolution. 

Our conclusions do not necessarily imply that the summation theorem is inherently 
false as an isolated mathematical proposition. Rather that it has been applied to the wrong 
context. This erroneous application arises due to two main disparities. In the first place, 
the argument for the summation theorem in the original Kacser and Burns paper (1973) 
is derived using a discrete formulation. However, in most subsequent work ( e.g. Kacser 
& Burns 1981 ) a continous formulation of the theorem is used. Hence the discrete proof 
for the summation theorem does not correspond to the continous form in which it is 
frequently used. In the second place, mutations necessarily involve discrete changes in 
enzyme activity. Hence, whether a continous version of the summation theorem holds or 
not ( which it would in cases where the flux function is homogenous, see Giersch 1988 ) 
does not say anything about whether a pathway will be robust or exhibit dominance with 
respect to discrete changes. In order for a pathway to be robust with respect to discrete 
changes, the discrete version of the summation theorem has to hold. In this paper we 
show that for dicrete changes of any size the discrete version does not hold, and hence 
that robustness is not an inherent property of a metabolic pathway. 

The matter addressed here has been a subject of debate for thirty years. In order to be 
clear about the objects that we treat and the mathematical properites that pertain to them a 
rigorous approach to the problem was necessary. The main assumptions and implications 
we derived are summarized as follows: 



INITIAL ASSUMPTION 

Assumption 0. The metabolic pathways studied are two-enzyme pathways. They have 
one input and one output and can reach a steady state flux. 



4 



CASE A : DIFFUSION LIMITED PATHWAYS 

Assumption 1. Maximal rate of input into the metabolic pathway is limited by a diffu- 
sion process. 

Implication 1. The flux summation theorem does not hold. 



CASE B : SEQUENTIAL PATHWAYS THAT ALLOW FOR ENZYME SATURATION 
Assumption 1. The two enzymes are near saturation. 
Implication 1. The flux summation theorem does not hold. 



CASE C : NUMERICAL ANALYSIS OF PATHWAYS WITH REVERSIBLE MICHAELIS- 
MENTEN KINETICS 

Assumption 1. Each enzyme in the pathway behaves according to reversible Michaelis- 
Menten enzyme kinetics. Isolated Michaelis-Menten enzymes do allow for enzyme 
saturation. Note that such kinetics are manifest in empirical studies of isolated 
enzymes in vitro. 

Assumption 2. When Michaelis-Menten enzymes are placed in the context of a multi- 
enzyme pathway, their kinetic mechanisms do not change. This is a pure assump- 
tion. 

Assumption 3. Maximal rate of input into the metabolic pathway is limited by a dif- 
fusion process. As in assumption 1, this is is an assumption that is empirically 
justified ( [Tralau et ai, 2000| ). 

Numerical Result 1. In regimes where both enzymes are saturated the summation theo- 
rem does not hold. 

Numerical Result 2. In regimes where flux is near the diffusion limited rate, the sum- 
mation theorem does not hold. 

Assumption 4. Mutations can cause heritable variation in the catalytic turnover rate k ca t 
of an enzyme. 
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Numerical Result 3. The control of a pathway can be switched from one enzyme to 
another via mutations affecting k ca t- Each enzyme can be the sole controlling 
enzyme of a pathway. 

The main reason for the restricted applicability of the summation theorem is the fact 
that for its construction the non-linear interaction effects on flux were ignored. It is pre- 
cisely the interaction effects between enzymes that lead to the manifestation of interaction 
effects ( epistasis ). Furthermore it is the interposition of the non-linear properties of en- 
zymes that allow for the evolution of control. 



2. Mathematical Introduction 

2.1. The conception of control in MCA 

Upon devising MCA, as Kacser and Burns put it, it was necessary to place the understand- 
ing of individual enzyme behavior into the context of a "metabolic society" composed of 



several enzymes (Kacser & Burns, 1979). In accordance to this goal, one of the objec- 
tives of metabolic control analysis was to ascertain what role each enzyme could play in 
affecting steady state flux in a particular metabolic pathway. Specifically, how important 
each enzyme is in its ability to affect steady state flux via the regulation of the enzyme 

concentration. Hence the use of the control coefficient as Cf = —/—=^, which is a 

non-dimensionalized sensitivity measure that shows the ratio of proportional change in 
flux to proportional change in enzyme concentration. 

2.2. The flux summation theorem 

One of the central tenets of metabolic control theory has been that the sum of the control 
coefficients in a pathway with n enzymes equals one. That is, 

n 

*£cf = l. (2) 

1=1 

Equation (||) is commonly referred to as the flux summation theorem and was derived by 
Kacser and Burns using the following form of reasoning: 

Given an unbranched chain of n enzyme catalyzed reactions, if we were to increase 
the concentration Ei of enzyme i by the proportional ratio 5Ei / Ei, the change in flux due 
to enzyme i would be given by rearranging (|l|) to give: 

For all enzymes i : ^ = (~]^J ' ^ 

If 5Ei /Ei = a , where a is a constant, one could consider the case where all enzyme 
concentrations are changed by the constant proportion a. In such a case Kacser and Burns 



6 



assume the following equation: 



5 jL - V (— 



(4) 



Hence, using (HI), (Eh and a = 5Ei/Ei ; 



5J_ 

T 



(5) 



i=l 



i=l 



Separately, the assumption is made that if all enzyme concentrations are uniformly changed 
by a , then 



SJ_ 
J 



a. 



(6) 



Substituting equation @ into (||) we obtain equation < 

n 



i=l 



For a clearer view, the Kacser and Burns derivation of (Q) can be summarized in the 
following conditionals: 



( (1); (Definition) [ (%) . = C % 
d); (Assumption) [ ^ = E?=i (¥) 



A 



A 



(7) 



] / 



\ ( Restriction ) (Vi) : [ 5Ei/Ei = a 



(8) 



(I); [^ = «E?=i^/] a \ 

( Restriction ) (Vi) : [ 5Ei/Ei = a ] 
V (D; ( Assumption ) [ ( (Vi) : [ 5^/£? 4 = a])^[f = a]]/ 

[^ = 1]). 



Note that technically even if the assumptions (|h and (||) were true, (Q) would be only 
proven for the restricted case where for every enzyme i : [ bEijEi = a]. 
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3. Construction of two-enzyme pathway models 



3.1. Two-enzyme pathways as dynamical systems 

A general dynamical model of multi-enzyme systems can be formulated based on the 
classical model of single-substrate enzyme catalysis. Under such a scheme the existence 
of an intermediate enzyme-substrate complex is posited and a set of forward and back- 
ward reaction rates is attributed to each transformation. Consider a pathway in which an 
outside substrate u diffuses into the initial reaction compartment,where the substrate is 
fed into a two reaction pathway that comprises two successive enzyme catalyzed reac- 
tions. An irreversible sink step is added to the end of the reaction sequence in order to 
remove the product. 



. , Diffusion Enzyme 1 Enzyme 2 Sink / , 

input(u) s\ <-> s 2 «-> S3 — > output(o) 

The above scenario is represented by the following kinetic model shown schematically: 



D 

u D Si 



fcl fc3 

ei + ai k2 (esi) k4 s 2 + &\ 



k5^ k7 t 

e2 + s 2 fc6 (es 2 ) fcs s 3 + e 2 



(9) 



s 3 — ► o 



where (es^) and (es 2 ) are tne enzyme-substrate complexes, while el and e2 are the first 
and second enzymes catalyzing the two respective reactions. The input and output cor- 
respond to u and o respectively. The kinetic model (^) is a class of model that can be 
represented as a system of autonomous differential equations with input u(t) and initial 
conditions s(0) and e(0) such that; 



s(t) = /i[s(t),e(t),u(t)] 

e(t) = / 2 [s(t),e(t)] (10) 
d(t) = Ms 3 (t)] 

where at time t, the vector s(t) = ( s\(t), s 2 (t), ss(t) } is the state vector of substrate con- 
centrations, the vector e(i) = ( ei (t), es^ (t), e 2 (t), es_ 2 (t) ) is the state vector of enzyme 
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and enzyme-substrate complex concentrations, the variable o(t) is the concentration as- 
sociated with the output variable and d(t) is the output flux. The notation x associated 
with any variable x denotes the time derivative dx/dt. If the system reaches a steady 
state flux the output flux 6(t) is equivalent to the steady state flux J used in the MCA 
literature. 

3.2. Characterization of two-enzyme pathways into distinct classes 

We will make the rest of our argument using a set theroetic framework. For any given set, 
if we can define the general properties of its members, then any object which is deemed to 
be a member of that set will also have those properties. A system of equations is a logical 
proposition. As such we can cast our dynamical system model in a set theoretic guise. A 
pathway p is an object. Each object has a series of properties K associated with it, such 
that K \p\ = True. The properties K delineate the context and system of equations that 
can describe a pathway ( For alternative representations see Fontana & Buss 1996 ). Each 
pathway p can then be treated as an element of a set with distinct properites. 

Definition 1 : Pi i = Set of two enzyme pathways with one input and one output. 

There are no strict stipulations on the stoichiometry nor on the specifics of enzyme kinetic 
mechanisms for membership in Pii ( see Appendix A ). 



Definition 2: P u ,j = Set of two-enzyme pathways with constant input and steady-state 
flux 

The set P u ,j is defined as the subset P u ,j C Pi i such that p € P u ,j if and only if the 
following four conditions hold: 

1. P€Pi,i 

2. Irrespective of the enzyme kinetic mechanisms, p has the stoichiometry of a two- 
enzyme sequential pathway as exhibited in model @. 

3. When driven by a constant input u, p exhibits a unique steady state output flux J* 
that can be associated with each pair of total enzyme concentrations (E\,E2). 

4. For any error margin e > 0, the steady state output flux J* can be approximated by 
using the total enzyme concentrations (E\, E2) and constant input it as arguments 
for a computable function g £ such that J = g £ [u, E\, E%\ and \ J — J*\ < e . 

Items 2 to 4 are elaborated more precisely in Definitions 2.2 to 2.4 in Appendix A. 
Note that the MCA approach holds informal versions of these propositions as premises, 
without which talking about the steady-state flux J and the summation theorem would 
not be possible. 
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Definition 3: P u d,j = Set of diffusion limited two-enzyme pathways with steady-state 
flux 

The set P u d,j is defined as the subset P u d,j C P u ,j such that p € P u d,j if and only if 
the following two conditions hold: 

1. P € P U ,J 

2. The net input flux from u(t) into the system is governed by a diffusion-limited pro- 
cess such that Netstartflux = D ( u(t) — si(t) ) where D is a diffusion constant. 

Item 2 is more precisely defined in definition 3.2 in Appendix A. 



Definition 4: P u ,v max ,J = Set of two-enzyme pathways with steady-state flux and kinet- 
ics that allows for enzyme saturation. 

The set P u ,v max ,j is defined as the subset P u ,v max ,j C P u ,j such that p e P u ,v max ,j 
if and only if the following two conditions hold: 

1. P G P U ,J 

2. For every enzyme i, there exists a constant k cat ^ > such that for all values of 
Ei > the steady state flux is consistent with the relation J < Ei k cat ^y 

Item 2 is more precisely defined in definition 4.2 in Appendix A. 



Definition 5: P u D,v max ,J = Set of diffusion limited two-enzyme pathways with steady- 
state flux and kinetics that allows for enzyme saturation. 

The set P u D,v max ,J is defined as the intersection 

PuD,V max ,J = (PuD,J n P u ,V max ,j) C P U; J. 



As an intuitive guide, the order of the generality of the representations discussed 
so far, in order of higher generality are: P u D,v max ,J C P u d,j C P u ,j C Pi,i and 
PuD,v max ,J C P u ,v max ,J C P u ,j C Pi,i (See Figure §. 
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4. Problems with the summation theroem 



The definitions from the previous section have a series of logical implications. Here we 
present these implications and explain their significance. The proofs for these proposi- 
tions are given in Appendix B. 

4.1. Pathways with a constant input and a steady-state flux (P u ,j) 
Proposition |0] 

For any pathway of class P u ,j. the flux summation theorem implies that for any given 
ratio of enzymes < (Ex/E^) < oo , there exists a Cg such that flux is a linear function 
of the magnitude of the vector (Ex, E 2 ). 

(Vp 6 P u ,j) : (V0 | tanfl = Ex/E 2 ) : (3C e ) : 

[ C( + C( = l J=\\(E 1 ,E 2 }\\C e }. {n) 



The meaning of proposition 11 is very simple. If a pathway were to be of class P u ,j 
and were to obey the flux summation theorem the following would hold: the flux would 
behave in a linear fashion with respect to tandem changes in which both enzyme concen- 
trations are changed by the same percentage. It is a short few steps (Appendix B) to show 
that 



V a 



SEx 5E 2 



Ei 



E 2 



y =a & J=\\(E 1 ,E 2 )\\C e 



(12) 



Proposition ( JT2h is equivalent to a statement of homogeneity ( See Giersch 1988 ), such 
that if J = g[Ei, E 2 ] then 



g[ (1 + a)E u (1 + a)E 2 } = (1 + a)g[ E lt E 2 



(13) 



For example consider two enzymes whose genes are located on the same operon. In 
addition consider the case in which such a cell is in a state in which the concentrations of 
the two enzymes are initially negligible. When the operon is turned on, each enzyme will 
have a particular rate of expression. If the ratio between those rates of expression remains 
constant, then the flux summation theorem would imply that the flux would increase 
linearly without reaching a saturation plateau. This is what J = Ce implies, 

where || (Ei, E 2 ) \\ = \JE\ + E\. Hence, in a pathway of class P u ,j, in which we have 
made no assertion about the particulars of the enzyme-kinetic equations or parameters, 
the flux summation theorem assumes a particular type of dynamic behavior. A flux that 
increases linearly without saturation is a very special case and such a condition can only 
hold in a restricted subset of P u ,j- In fact, as proposition 14 will attest, even for a 
simple pathway with the sole kinetic imposition of a diffusion barrier, the flux summation 
theorem fails. 
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4.2. Pathways with a diffusion barrier (P u d,j) 
Proposition |l4| 

For any pathway of class P u d,j» the flux summation theorem is false. 

(Vp G P U D,J) : [ (V<5£ x , <5£ 2 ) : {^E U E 2 G M^ ) : [ Cf + C 2 J ^ 1 ] ] (14) 



Proposition 14 basically states the general result that any pathway which has a diffusion 
barrier at a point before its first step will not obey the summation theorem. This includes 
almost any metabolic pathway that takes place in the confines of a cell. In fact it can be 
shown that any pathway that exhibits a plateauing effect on its flux surface will not obey 
the summation theorem. The proof of this theorem is intuitively simple (Appendix B). It 
is basically dependent on the conditional relation 

C( + Ci = l => J=\\(E 1 ,E 2 )\\C e (15) 



from (11). For any pathway class that is a subset of P u ,j, if the pathways do not obey the 



right hand side of [15| then they do not obey the left hand side either. For P u d,j one can 
show that J = || (Ei,E 2 ) \\ Cq does not hold when Cq > 0. This is because the diffusion 
barrier will impede flux from increasing indefinitley. This implies that Ya=i = ^ 
cannot be true for all enzyme concentrations. This falsehood applies to all diffusion- 
limited flux surfaces except for a disfunctional cell with flux surface J = 0. Proposition 



14 is due to the fact that control is not only dependent on the enzymatic steps in a pathway, 



but also the diffusion steps. 



4.3. Pathways that allow for enzyme saturation (P u ,v raax ,j) 
Proposition |6| 

For any pathway of class P u v rriax .J. if both enzymes approach saturation and changes in 
enzyme concentration are positive, then the sum of control coefficients approaches zero. 



(Vp G P u ,v max ,j) 



(V5E U 5E 2 G R^o) : lim {C( + C J 2 ) = 



(16) 



Whenever two enzymes are near saturation, a situation is created whereupon the flux ef- 
fects of increasing the concentration for one enzyme is restrained by the other enzyme. 
Hence whenever two enzymes approach saturation simultaneously, the sum of their con- 
trol coefficients approaches zero and the flux summation theorem does not hold ( see 
appendix B ). In fact it is not too difficult to extend the proof such that for a sequential 
pathway of any length, the sum of all control coefficients will approach zero whenever 
any two enzymes approach saturation. 
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Proposition |n| 



For any pathway of class P u ,v max ,J in which a decrease in enzyme concentration in- 
creases saturation, if both enzymes approach saturation and changes in enzyme concen- 
tration are negative, then the sum of control coefficients approaches two. 



For all intervals a\ < Ei < bf. 

If P G P u ,v max ,J and Mi : ( > and a { > and SEi < ) then: 

« 1 1 SU 1 ' C 'i' + ^)= 2 - <17) 



In such situations the pathway is no longer robust to decreases in enzyme concentration. 
Hence dominance is not an inherent property of these pathways. It is simple to extend 
this result to sequential pathways of any length, in which the sum of control coefficients 
will approach the the number of enzymes which are approaching saturation. 

r 



From propositions 14 



false in pathways of class P u D,v 



and 16| it also follows that the summation theorem is also 
J- 



4.5. Enzyme interactions in P u J 

The motivation for analyzing multi-enzyme pathways is that they behave differently from 
single-enzymes. The natural question is to ask what are the direct or indirect interactions 
between enzymes that lead to the difference between single-enzyme and multi-enzyme 
pathway behavior. Interaction does not necessarily have to signify direct physical interac- 
tion between two enzymes. It can occur via the effect of one enzyme on the concentration 
of pathway intermediates that in turn physically interact with another enzyme. In equa- 
tion (6) which was one of the assumptions that Kacser and Burns made in deriving the 
summation theorem, the proposition ^j- = Ya=i (t') • amount s to an assumption that 
enzymes do not exhibit interactions. Here interaction is specifically referring to the mu- 
tual effect that the two enzymes would have on each-other's ability to affect pathway 
flux: 



Proposition |^ 

For any pathway of class P u ,j, the proposition [ = Y^i=i irr) ■ ] ® an ^ ^ e re ~ 



striction a 



SEi 
Ei 



SE 2 
E 2 



imply that the change in flux due to change of total enzyme 
concentration at one locus is independent of total enzyme concentration at the second 
locus in the pathway. 



(VpeP u ,j) : (V£i,£ 2 ERl ) 



<9 2 J 
dE x dE 2 



(18) 



The assumption of independence of effects is something that surreptitiously goes against 
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the motivation for developing MCA. Nonetheless the assumption of no enzyme interac- 
tions is built into the summation theorem from the beginning. In fact: 



Proposition |l^ 

For any pathway of class P u j, the summation theorem implies that the change in flux 
due to change of total enzyme concentration at one locus is independent of total enzyme 
concentration at the second locus in the pathway. 



(Vj9GP u ,j) : (V^i.Ea G 



[(VS£i, 8E 2 ) ■ C{ + Ci = l 



d 2 J 
dE-i dE 2 



(19) 



This means that the summation theorem implies the inexistence of epistasis. It immedi- 
ately follows that 

(Vp € P uD ,j) : V (E X , E 2 G Ri Q | + ) : 

[ (35E 1 ,SE 2 ): [Cl + Ci^l] ]. 

Hence in any regime where epistasis is present the summation theorem is violated. Fur- 
thermore: 

Proposition [2l] 

For any pathway of class P U ,J. the summation theorem holds if and only if flux is a linear 
function of enzyme concentrations. 

(VpGPu ( j):(V£i,i? 2 eRg: 

| [(VSEi, SE 2 ) : C( + C$ = l] J = E x c x + E 2 c 2 ] . V ; 



In short, the summation theorem imposes conditions on enzyme interactions that are 
neither explicitly argued for nor conceptually justified. 

5. Numerical analysis of a pathway of class P U D,v max ,j 

Here we use numerical simulations for a particular case in order to verify the analytical 
results from the previous section. Consider a case where the enzymes exhibit reversible 
Michaelis-Menten kinetics. A diffusion barrier is included at the beginning of the path- 
way. During the interval of analysis < t < tf, the input remains constant such that 
u(t) = (e.g. approximation of large nutrient pool and/or steady state conditions of 



a chemostat). In such a case the full form of the system of differential equations Q10|) 
combined with the kinetic model (0) are: 



14 



' dt 

d\eS] 



-D [ai] - fcl [ei] [si] +D[u}+ k2 [es x ] 



=li = A;l ei [si] - A:2 [es.i] - k3 [esj + M [ei] [s 2 ] 



d[s 2 ] 



^f = k3 [es.i] - £4 [ei] [s 2 ] - A;5 [e 2 ] [s 2 ] + A;6 [es 2 ] 



dies 



=jr i = k5 [e 2 ] [s 2 ] - k6 [es 2 ] - k7 [es 2 ] + k8 [e 2 ] [s 3 ] 



d[s 3 ] 



V if = kl fe&l - Q M - k8 [e 2 ] [s 3 ] 

-fcl [ei] [si] + A;2 [es x ] + k3 [es x ] - &4 [ei] [s 2 ] 



d[ei] 



d[e 2 ] 
d[o] 



-fe5 [e 2 ] [s 2 ] + A;6 [es 2 ] + kl [es 2 ] - k8 [e 2 ] [S3] 



(22) 



/ 



where the postfix (t) for each state variable is omitted. 

The system of differential equations ( |22| ) will serve as the basis for the examination 
of a particular case. For the present analysis we are interested in the steady states of the 
dynamical system. At steady state, as defined in appendix A, the derivatives s(t) and e(i) 



are equal to zero. Consider a case where the equations in ( J22| ) serve as the basis for the 
functions fx, f 2 and / 3 with the stoichiometric constraint that 



[Ei] = [ex] + [esx] 
[E 2 ] = M + [es 2 ] 



(23) 



In such a case, the steady state solution can be obtained in an implicit form ( See appendix 
C ). The implicit equations can then be solved numerically by Newton-Raphson based 
algorithms that converge to the flux surface shown in Figure ^[ On a cursory level the 
surface shows the plateau effect familiar in most models of two-enzyme pathways. The 
plateau arises due to the limiting effects of the diffusion barrier. We will refer to the 



pathway based on equations ( |22| ) as pathway pj^J. It is not difficult to discern that pathway 
is a member of the set P u D,v raax ,J- It satisfies the conditions required in definitions 
1 to 5. 



5.1. Test of the summation theorem 

Given the flux surface on £ P u D,v max ,J. the flux summation theorem can be tested 
directly. Figures ||] and |] show surfaces for for discrete versions of the sum Ya=1 Q 7 f° r 
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It is evident from these examples that in the pathway p|| there is no evidence of 
an invariant relation in which the control coefficients sum to the predetermined value of 
one. Figures | and | serve as illustrations of propositions |l4|, [l6| and [n|. The summation 
theorem fails in the regimes where flux is near its diffusion-limited maximal rate or in 
areas where both enzymes are near saturation. The regions where both enzymes approach 
saturation are shown in Figure ||| 



5.2. Epistatic interactions between enzymes 

The question of the extent to which enzymes have interactive effects on flux can also be 
addressed directly using the pathway p|| as an example. In Figure |6| the mixed partial 
difference 5 2 J / 5E\5E2 is used as a measure of epistatic interactions between enzymes 
(proposition [18|). Note that the region where the epistatic interactions increase coincides 
with a region where the summation theorem fails and both enzymes are approaching 
saturation. 



5.3. Rate limiting enzymes and the switching of control between enzymes 

One of the claims from MCA has been that there is no such thing as a rate limiting 
enzyme in a pathway. The associated claim is that control is mostly shared between 



enzymes in a pathway ( |Kacser, 1995[ ). The results of our numerical simulations of ; 
do not support this claim. Figure [7] shows the numerical results for two pathways with 
different k cat configurations. In the range of concentration regimes {E\, E2) examined, 
most regimes manifest a sole enzyme as the controlling enzyme for a given pathway. In 
regimes where flux is near the diffusion limited maximal flux, no enzyme has control. The 
only region where the enzymes share control is where the summation theorem fails. In 
fact in the regions that the summation theorem holds, the sum Ym=\ = 1 i s basically 
the value C% for the sole controlling enzyme z at that regime. Furthermore the identity of 
the controlling enzyme at a given regime is not invariant. The identity of the controlling 
enzyme can change by changing the k cat values of the enzymes involved (e.g. by mutation 
). In Figure ^ the switch in control is clearly apparent by comparing the values of C% 
between the two k cat configurations. 



6. Discussion and conclusion 

The main intent of the present paper has been to present an analysis and critique of some 
of the main biological conclusions derived from MCA. Our results strongly support ini- 
tial objections that were made with regards to the MCA perspective on the biochemical 
nature of dominance ( Cornish-Bowden 1987, Savageau 1992. See also Mayo & Burger 
1997 ). At the core of the problem of multienzyme systems is interaction effects between 
enzymes. In order to address these interactions, the non-linear properties of enzyme catal- 
ysis cannot be ignored. These non-linear interactions play a major role in the phenotypic 
characteristics of metabolic physiology. Our analysis shows that for discrete changes of 
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any magnitude, the flux summation theorem does not hold. This means that there are no 
a-priori constraints that would require the magnitude of mutational effects to be in a low 
range. This leaves open the possibility for two phenomena which had been previously 
rejected in MCA. In the first place, rate limiting steps or controlling enzymes do not have 
to be rare, and their identity can change by evolution. Secondly, phenotypic robustness 
with respect to mutations is not an automatic property of metabolic pathways; the de- 
gree of robustness can be modified through the evolution of kinetic properties. These are 
implications that we will examine in a later work. 
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Appendices 
A. Definitions 



In this section some of the conditions proposed in definitions 1,2, 3 and 4 are formu- 
lated in a more precise language as definitions 2.2, 2.3, 2.4, 3.2 and 4.2. In these def- 
initions quantifiers are encapsulated in parenthesis and propositions in square brackets. 
The symbol "A" denotes the boolean expression "and". The symbol "M> " denotes the 
real numbers greater than or equal to zero. Ordered tuples and vectors are encapsu- 
lated by the symbols "( )". The property K associated with an object p is denoted as 
("name of property") [p\ . 

Definition 1 



(Vp € Pi,i) : p 



s(t) = h[s(t),e(t),u(t)] A 
e(t) =/ 2 [s(t),e(t)] A 
6(t) = f 3 [s 3 (t)] A 
(u(t),o(t) eRl Q ) A (e(0) G 



^ 4 >o) 



A (s(0)€R| o ) 



(24) 



Definition 2.2 

The proposition ("Stoichiometry: two enz. seq.")[pj imposes both organizational and 
kinetic constraints on the pathway p G P u j by designating the stoichiometric proper- 
ties of the pathway. Such a proposition is necessary for an abstract representation of 
physical mass balance. 



("Stoichiometry: two enz. seq.")[pj = 



(V(t,u(t),o(t)) Gig: (Ve(0) G K 4 >0 ) : (Vs(0) G 



V v 



[E 1 (t) = E 2 (t)=0] A 
[E 1 = E 1 (t) = e 1 (t) + es 1 (t) 
[E 2 = E 2 {t)=e 2 (t) + es 2 (t) 
[E 1 >0 A E 2 > ] 



A 
A 



A 



f(t) = ai(t) + 2es 1 {t) + s 2 (t) + 2es 2 (t) 
+ a 3 (t) +ei(t) + e 2 (t) 



[ e(t) = / in (t) - Wt) - 6(t) 



v >0 



A 



(25) 



where the pair (Ei,E 2 ) G M> are the total enzyme concentrations for any pathway p 
such that (E\,E 2 ) does not change for each particular vector of starting conditions and 
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that for a given set of initial conditions, the concentrations E\ and E 2 are to be deter- 
mined at any time t by summing the free-enzyme and enzyme-substrate complex concen- 
trations for each enzyme respectively. The variable £(t) needs some explanation. £(t) 
determines the "stoichiometric content" of the pathway in the sense that the coefficients 
Cj in the equation 

£(t) = cisi(t) + c 2 es 1 (t) + c 3 s 2 (t) + c 4 es 2 (i) + + c 6 e x (i) + c 7 e 2 (t) 

are determined by whether in reference to the input molecule, the transformations leading 
to the macro-molecules associated with each coefficient are fusions, cleavages or one-one 
transformations. The values of the c, in (25) correspond to the kinetic model (0). The 
equation for also deserves attention. Each metabolic pathway can be characterized 
such that £(t) is determined by the additive combination of two values at time t: 

= Netstartflux(t) - Netendflux(t) 



which in turn can be subdivided as 



Netstartflux(t) = h n {t) 
Netendflux(t) =0 ou t(t) 



O rev (£) 



and 



where Ii n {t) is the flux entering the system from the input variable u(t), I rev (t) is the 
flux exiting the system and going back into u(t), ou t(t) is the flux exiting the system 
into the output variable o(t) and O rev (t) is the flux entering back into the system from 
o(t). For a kinetic model with an irreversible sink step, O rev (t) = 0. For a system 
where the output variable o(t) is only driven by a single pathway p, ou t(t) = d(t). 
Finally for P u ,j the only stipulations made for the two components of Netstartflux(t) 
is that they must both be greater or equal to zero such that Ii n (t), I rev (t) G M> and 
Netstartflux(t) = Ii n {t) — I rev (t). 



Definition 2.3 

For a pathway p G P U ,J the logical proposition 

("Const, u — > Unique (u, E\, E 2 , J*)") \p\ states that given an input u(t), if u(t) re- 
mains constant then the pathway will reach a unique steady state flux J* G M> that 
is dependent on total enzyme concentrations E\ and E 2 such that: 



Unique (u,Ei,E 2 ,J*)")[p\ 



("Const, u 

(Vn(0) £ 1R>q ) : (V^i(0), E 2 (0) G M^ ) : 
(3J* eRl ) : (VoWeRio) : (V Sl (0), s 2 (0), a 3 (0) G 

(0,0,0) A 

(0,0,0,0) A 

A 

J* 



(Vt Gl 
[«(*) 



'>o) : 



V 



lim s(t) 



lim e(t) 



lim d(t) 



lim d(t) 



(26) 



Note that at steady state £ (t) = 0. 
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Definition 2.4 



The proposition ("computable: (E\,E 2 , u) i— > J ") [p\ states that an approximation J of 
J* is computable as a function of the initial input and total enzyme concentrations: 

( "computable: (Ex ,E 2 ,u) <—>■ J")|_pJ = 

(V(u(0),o(t)) € 4 ) : (Ve(0) G R| ) : (Vs(0) G R| ) : 
( [J = g £ [u(0),E 1 ,E 2 ] A |J- J* | <e] ) 

The above definition means that for any error margin e > 0, there exists at-least one 
effectively computable function g £ that maps to the approximate value J by using the 
total enzyme concentrations and initial input as arguments. Note that g £ as an effectively 
computable function refers to both numerical and analytic functions that can produce an 
approximation J. The assumption of computability is associated with the assumption 
that there is a rule based mechanistic relation between E\ , E 2 and u. There is a strong 
assumption inherent in the definitions of g £ in P U) j: that for any pathway p G P U ,J a 
steady state flux exists and can be determined using only u(0), E\ and E 2 as arguments. 
That is, except for the constraints imposed by n(0), E\ and E 2 the steady state is indepen- 
dent of the specifics of starting conditions. Note that for classes of pathways that exhibit 
multiple equilibria the trajectories become important and all the starting conditions have 
to be included as arguments in a newly defined computable function other than g £ , which 
would be embedded in a system other than P u ,j- 

Definition 3.2 

( "Diffusion Limited 7" ) [p\ = 

(3D G Ri ) : 

(V(t,u(t),o(t)) G M> ) : (Ve(0) G M| ) : (Vs(0) G R| ) : (28) 
/ [I in (t)=Du(t)\ A \ 

V [I rev (t) = D Sl (t)\ J 

where D is the diffusion constant. 
Definition 4.2 

("Saturable Enzymes") [p\ = 

(3(k cat(1) ,k cat{2) ) G Ri ) : (V(E U E 2 ) G Mi, ) : 
(V(t,u(t),o(t)) G^o): (Ve(0) G R| ) : (Vs(0) G R| ) : (29) 
/ [ J < ^iA; cat(1) ] A 

V [ J < E 2 k cat{2) } 
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B. Proofs 



B.l Proof of proposition [Tl| 

We start with the summation theorem Y^7=i = Using the definition of Cf from ([I]) 
we have 



(yE u E 2 e 



5 J Ex 5 J E 2 _ i 



5^1 J <5£ 2 J 



The continuous version of (|^) is ^ lim ( £i + j^- E 2 



5Ei,SE2-*0 



(30) 

J which is equiva- 

dEx " x ' 5^2 " z (31) 
Equation is a first order quasi-linear partial differential equation with general solution 



lent to 



9J „ dJ „ 

Ex + — — E 2 = J. 



J = E 2 f 



Ex 
E 2 



Ei f 



El 
Ei 



(32) 



where /, /' : M 1 i — > R l (see for example Rade & Westergren 1995, p230). For each fixed 
ratio of Ei/E 2 , the vector (Ex, E 2 ) forms a fixed angle 9 with the coordinate axes E 2 
and Ei such that E%/E 2 = tan 9. Hence we can formulate the following proposition: 



(3g st : M| » R 1 ) : (VEx,E 2 G K^) : (3f st : E^ h- R 1 ) 
J = g st [E 1 ,E 2 ] = E 2 f st [t&n9] = h cos 9 f st [tan 9] 



where h is the norm of the vector (Ex,E 2 ) such that h = \\ (Ei, E 2 ) \\ 
Given that for every 9 there exists a Cq such that Cq = cos 9 ■ f s t [tan 8] , we have 

(V0 \ 0<9< it/2) : (3C e ) :[J = hC e }. 



(33) 



VW+ei. 



(34) 



proposition |11 



B.2 Proof of Proposition |g 

The proof of the proposition 



a 44> 



J = ||(^i,^2)|| Cq] is uncomplicated. We 



start with proving the forward conditional [" = a => J = [| (£1 , E 2 ) \\ C e ] . By defi- 



nition we are given that 



a = 5 Ex /Ex = 5E 2 /E 2 . 



Consider the vector (Ex,E 2 ) and its norm h = \\(Ex,E 2 



(35) 

^E\ + E\ . Given 



simultaneous increases in enzyme concentration 5 Ex and 5E 2 that obey (35) we have 
h + 5h = + a) 2 ^ 2 + (1 + a) 2 El. Hence = a and thereby 



<J J/ J = Sh/h. 



(36) 
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In its continuous form equation (|g) can be set up as a separable differential equation 
such that 



Integrating both sides of (|37|) we have the solution 



dJ dh 
~J = ~h 



(VE 1 /E 2 ) : (3Cg) : [J = hCg] 



SJ 



The backward conditional [ 
entiating (^8|) with respect to h. 



a 



(37) 



(38) 



J = || {Ex, E 2 ) || Cg] can be proven by differ- 
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B.3 Proof of proposition [1^ 

Consider a case where we posit the existence of a pathway px G P U .J which obeys the 



flux summation theorem . In such a case, using (34) we have; 



(V0 | < 9 < it/2) : 
[[h—>oo] o [ Ex —> oo A E 2 



oo 



(39) 



Furthermore for the function g st defined in ( |33[ ) we have the following implication: 

{{3 P1 G P UiJ ) : [(VEx,E 2 G R l >0 ) : £™ =1 Cf = l] ) : (V£ x , £ 2 G M^ ) : 
lim J = lim g s t[Ex,E2] = lim /i • cosci • / s t[tanci] = lim h Cg = oo. 

h— >oo h^oo h^oo h^oo 



(40) 



On the other hand consider a diffusion-limited pathway p 2 6 P u d,j where J = , £"2] 
as defined by (^l|). When the system reaches a steady state as defined by (f26|), we know 
by (EBh that £(t) = and = - I re „(i) - d(t). Since lim 6(f) = J*, using © 

t^oo — 

we have the steady state condition: 



Since u(t) = 0; 



J* = D u(t) — D sx(t) and 
J±e = Du{t) - Dsx{t). 



(Vp 2 G P uD ,j) : (31?) : (Vn(0) = c 2 ) : (Ve) : 
[c 3 = Dc 2 ] A [ J ± e = c 3 - Dsx(t) ]. 



(41) 



Recall that D > by (g8|) and > by @. For any finite < e < 00 we can 
define a constant C4 such that 04 = 03 + e. Consequently 



J < C4 



(42) 
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where we know that < C4 < 00. Hence 



(yP2 ePuD,j) : (V(£i,£ 2 ) eKj> : 
lim J = lim g[E\,E^[ < (04 < 00). 

h— »oo h—>oo 



(43) 



If we now posit a pathway y»3 that both obeys the flux summation theorem and is diffusion- 
limited, propositions (fw|) and (43) lead to the proposition: 



3(p 3 e P u d,j) : 



P3 G Pu,J ] A [(yEk,E2 e 4 ) : E?=i C/ = 1] A 

[ OO < (C4 < 00) ] 



(44) 



which is a contradiction. Therefore 



(Vp G PuD,j) : 



i=l 
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B.4 Proof of proposition [R| 



For a pathway p G Pu,v raax ,J, consider a case in which enzyme 2 is saturable such that 

(45) 

where for any 6 Ex , J' = J + S J. It follows that if E2 is held constant, for any SEx 



(V^iGM^o): [J<E 2 k cat{2) A J' < E 2 k cat{2) ] , 



5 J E 2 k cat{2 ^ - J 



8E1 ~ SEx 
Multiplying both sides by E\ / J we have 



E x SJ EiE 2 k cat{2) - J Eh. 



J SEx 



JSEx 



Hence by ([[]) we obtain 



j ExE2k cat (2) Ex 



a < 



(46) 



(47) 



(48) 



JSEx SEx ' 

In a similar fashion if enzyme 1 is saturable and Ex held constant, for any SE 2 we have 



j , ExE 2 k cat (x) E 2 
JSE 2 ~SE 2 



Ci < 



For each enzyme i, consider a measure of saturation Sati such that 

Sati = J / Eik cat ^y 



(49) 



(50) 
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It follows that for any enzyme i, 



lim 

5otj-+l 



E,k 



i K cat(i) 
~J 



1. 



(51) 



Denoting R.H.S.($tf) as the right hand side of (48), as enzyme 2 approaches saturation, 
@ and @ imply that 



sii ( M . s . @) = J.-^ = o. 



Similarly, if enzyme 1 is approaching saturation, ( |51| ) and ( |49[ ) imply that 

i? 2 -£?2 



(Vp G Pu,V raax ,j) 



lim 

Sati ,Sat2- 



(52) 



(53) 



Consequently when both enzyme 1 and enzyme 2 approach saturation, (p2|) and (|53|) 
imply that 



(54) 
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B.4 Proof of proposition [17] 

For a pathway p G P u ,v raax ,j> consider a case in which enzyme 2 is saturable such that 



(V^i€M>o): [J<E 2 k cat{2) A / < fyk^m ] , 

where for any 5E 2 , J' = J + 5 J. 

Given that SaU = J/Eik cat ^, if E\ is held constant, for any 5E 2 < then 

lim f E2k ™t( 2 A = i and lim (J) = £2^(2) • 



(55) 



(56) 



Given (^) and (|56|), for domains in which decreasing E 2 increases saturation we have: 

For all intervals a < E 2 < b : 

If > and a > and 8E 2 < then: 



lim iS J) = 5E 2 k cat{2) . 



Hence 



lim ( C j) = Mm ( = s I^ml^m = 1 

at 2 ^l V 2 > Sah-+1\6E/Ej 5E 2 /E 2 



Sat 2 



(57) 



(58) 
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Similarly if enzyme 1 is saturable we can show that 

lim (C() = 1. (59) 

Hence for all intervals a» < E; < 6, : 



(Vp G Pu,V max ,j) : 

(60) 

lim (C( + C() = 2 
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B.5 Proof of proposition ^ 



The most straight forward proof of (18) is by direct deduction. To do this we define a 
set of notations for representing discrete changes in steady state flux due to changes in 
enzyme concentration. 

Definition B.4a: 

For a given input u(t) that is constant, let J be the abbreviated notation for the result of a 
mapping in a function g : R 2 i— > R 1 , with the variables E\ and E% as arguments: 

(Vp G P u ,j) : 

(V«(0) G R 1 ) : (V(e > 0) G R^ ) : (3g : R| ^ R 1 ) : (61) 
[J = g[E 1 ,E 2 ] & J = g £ [u(0),E 1 ,E 2 ]} 

The discrete operator ( ) applied to J is defined as: 

jf^?= g[Ek + SEi , E 2 + 5E 2 ] 

where the symbol "=" denotes identity by definition and for a given function g, the incre- 
ments in E\ and E 2 are held constant such that 5E\ = c\, 5E 2 = c 2 and c\,c 2 G R> - 
Similarly, 

Jh^= g[Ex +8E X , E 2 ) and J ^> = g[Ei , E 2 + 5E 2 }. 

Note that the operator ( ) i— can be applied iteratively and is also commutative in its 
iteration order such that 



25 



Definition B.5b: 



Let the operator ( )(V 



() 



applied to J ^— ? be defined as 



J 



E 1 E 2 . 



- J 



where from Definition Al we get (J 
Similarly, 



-J) = + 5£ 1 ,£ 2 + £E 2 ] -5^1,^2]). 



-J and J i — >= J i — * 



-J. 



Proof Argument: 

To prove the relation 



°>o) 



(V£7i,£ 2 G 



d 2 J 
dEx 8E 2 



0. 



(62) 



we start with the left hand side of the relation, which is equation (Q) from the Kacser & 
Burns derivation. Using definitions B.5a and B.5b for a two enzyme pathway, equation 
(H) is equivalent to 



J 



EiE 2 



J 



Ei 



J 



E 2 



J 



J 



+ 



J 



(63) 



By cancelling the denominator in (|63|) we proceed with a series of deductions in which 
the horizontal symbol denotes the biconditional relation " if and only if " between 
successive lines. The explanation for each biconditional is given to its right: 
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(VE 1 ,E 2 eRl ) : Vf a 



SEx SE 2 



E x 



J 



E 2 

E 1 E 2 



A A 



$ using Definitions B.5a & B.5b 

g[Ex + SEh , E 2 + 5E 2 ] - g[Ex,E 2 ] = {g[Ex + 5Ex , E 2 ] - g[Ex,E 2 ]) 

+ (g[Ex , E 2 + 5E 2 }-g[Ex,E 2 }) 
$ rearrange 

(g[Ex + SEx , E 2 + 5E 2 ] - g[Ex + SEx , E 2 \) - (g[Ex , E 2 + 5E 2 ] - g[E x ,E 2 \) = 

$ using Definitions B.5a & B.5b 



(g[Ex , E 2 + 5E 2 ]-g[Ex,E 2 }) 



(VEx,E 2 G 



^>o) 



V a 



SEx _ SE 2 
Ex E 2 







$ using Definitions B.5a & B.5b 



0. 



Hence 



SE 2 
E 2 



{WE 1: E 2 GRio) : V(a = ^ 

E 2 s E x 







(64) 



A A 

is operationally equivalent to the sentence in propo- 



The proposition ((J) 

sition [l4] that proposes the independence of control effects. Next it is simple to show 
equivalence to a continuous version such that 



(VEx,E 2 €Rl ) : 



V(a = ^ = ^): ((J) 



SE 2 
E 2 



E 2 







d 2 J 
Ml dE 2 







(65) 



From the third line of our sequence of biconditional relations we know that the left hand 
side of (|65|) has to comply with the equivalence 



(66) 



(V^^GRio): V(a = fi = ^j : 

((J)3V)h^=0 o -J^-J^+J = 

vv ' A A 

We also know by the general definition of partial differentiation that 



dJ 



lim 



J 



-J 



dEx SEx^o 



5Ex 



and 
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Substituting (56) into ( |67| ) we get (65). Using (65) and (|64|) we get ( |62|). 

■ proposition |l8| 



B.6 Proof of proposition [1^ 

From (11) and ( |l2| ) we know that the summation theorem implies that 

J 5E 1 5E 2 \ 5J h2 
V la = — — = — — : — -r— = a. 



Ei Ei 



J 



Let 



6E\ = aE\ and 5E 2 = aE 2 . 



By definition we know that 



{VE U E 2 G Ri ) : tySEi, 5E 2 G R^) : 

Ei=i c '/ = 1 J 



Using (68) and (70) we have 



c j T SJ\ 5J 2 

dJi2 = aJ = ——aE 1 + ——aE 2 . 

6E\ 6E 2 



Using (71) and (69) we have 



<JJ 12 = 5J 1+ SJ 2 . 



By (65) in the proof for proposition 18 we know that 



(VE 1 ,E 2 £R 1 >0 ) : 



v < a = it = it ) : 6 J ^ =SJi + SJ 2 



d 2 J 

dE-i dE 2 
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B.6 Proof of proposition 121 



Let J = J[ E\, E 2 ]. Given proposition 11 we know that 

J = hC e and J[0,0] = 0. 
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By@: 



Let 



V^iGR^: J[E 1 ,0]=E 1 C B2= o. 



5J ^Ei[E 1 ,E 2 } = J[E l ,E 2 ]-J[0,E 2 }. 



Using (|74j) and ( |75| ) we have 



VEi G rol • 



<yjo- B1 [^i,o] = j[£?i,o] - j[o,o] = ^iC E2=0 . 



Separately proposition [19] implies that 

6Jo->ei[Ei,E 2 ] = c A 



\/Ei G Rl Q : 3c A : V£ 2 € M> : 



Using (76) and (77) we have 



(VEi,^ GR| ) : 

^0->£l[ ^1) ^2 ] = E\Ce2=0- 



Separately using a procedure similar to deriving ( |75[ ) we have 

V£ 2 GKi : ^[0,£ 2 ] =E 2 C E1=0 . 

Given 

J[^i,£? 2 ] = J[0,£ 2 ] + 5J ^i[£a,£ 2 ] 



and equations (|78|) and ( |79| ) we are left with 

J[ Ei, E 2 ] = E 2 Cei=o + E\Ce2=q- 



(75) 



(76) 



(77) 



(78) 



(79) 



(80) 
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C. Derivations for Numerical Analysis 
Derivations for flux surfaces 

At steady state the algebraic solution for the state variables of the equations in (^2|) can 
be obtained using algebraic solving routines. The solution for the seven internal state 
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variables are in the form of five equations in which there are no explicit solutions for e± 
and e 2 such that: 

si = (e2Q(k2 + k 3 )k 5 k 7 + eik 2 k 4 (Q(kQ + k 7 ) + e 2 k G k s )) @ 

Mi = eiki(e 2 Qk 5 k 7 + eik 4 (Q(k 6 + fe 7 ) + e 2 k e ks)) Q 

s 2 = e 1 k 1 k 3 {Q(k6 + k 7 ) + e 2 kQk 8 )@ (81) 

es 2 = eie 2 kik 3 k 5 (Q + e 2 /cs) © 

«3 = ei e 2 kik 3 k 5 k 7 @ 



where 
G = 



uD 



e 2 Q{k 2 + k 3 )k 5 k 7 D + e\{Qk 2 k 4 {k & + /c 7 )D + e 2 (Qk\k 3 k^k 7 + k 2 k^k 6 ksD)) ' 

(82) 



In addition recall the stoichiometric constraint that 

(83) 



Ei = ei + esi 



#2 = e 2 + es 2 - 



Taking the equations (81), (|82|) and ( |83[ ) we have a system of seven equations. These 
can be solved by a convergent Newton-Raphson based numerical algorithm for any case 
where the set of real-valued parameters 

{D, E 1 ,E 2 , Q, h,k 2 , k 3 , fc 4 , k 5 , k e , k 7 , k 8 } 

is given. The solution for {si, s 2 ,s 3 , e\ , es 1 , e 2: es_ 2 } that is arrived at by numerical ap- 
proximation can then be used to solve for the steady state flux J. The surfaces in Figure 
U and Figure ^ are numerical solutions of (81 ), (|2|) and (p3|). 



The wild-type kinetic parameters used in the numerical procedures were k\ = k^ = 
4 x 10 7 M" 1 sec~ 1 ,k 2 = k 6 = 4x 10 2 sec~ 1 , k 3 = k 7 = 7 x 10 2 sec _1 and k 4 = k 8 = 
1 x 10 6 M _1 sec _1 . The diffusion constants used are D = Q = 3 x 10 1 sec~ 1 . For each 
enzyme K eq = 70 and for the wild type k cat = 42, OOOmm -1 . Figures || and ^ show the 
steady state solution of all state variables in the pathway dH. 

For simulations of mutations affecting k cat consider enzyme 1 as an example. Kinetic 
mutations have to occur under the thermodynamic constraint that the equilibrium constant 
K e q = k\k 3 /k 2 k\ for the reaction cannot change. Note that V max = k cat [Ei] = k 3 [E±\. 
For instance a reduction of k ca t to 25% of wild type can be modelled by multiplying each 
of the wild-type rate constants k 3 and k± by 1/4. 
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Failure of a modified summation conjecture 



A possible counter to the failure of the summation theorem is to try to include all steps in 
a pathway, including the non-enzymatic diffusion steps. Hence a conjecture that 

n 

J2cf = C J D + C J El + C J E1 + C J Q (84) 
i=i 

Here we show the numerical evidence that such a conjecture does not hold in In any 
case if such a conjecture had been true, it would not have been of much practical value 
since it would mean that for any pathway deep in metabolism one would have had to 
take into account every step starting from the initial diffusion barrier. A modification of 



the proof for proposition ( |16[ ) can be used to show the failure of the modified summation 
conjecture in regions of saturation. In Figure [K] we show the numerical failure of the 
modified summation conjecture for a case where 5E{ = 0.02^M and 5D/D = 5Q/Q = 
0.01. 
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PuD,Vinax.J 

Figure 1 : The relation between different classes of two enzyme pathways. 
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[ EI ] 




Figure 2: Metabolic flux as a function of total enzyme concentrations. The constant en- 
vironmental input is u = 0.75 mM. The catalytic turnover rates are symmetric such that 
k ca t(i) — k ca t(2) — 42, OOOrain -1 . Enzyme concentrations [Ei\ in fiM. Flux in mM sec^ 1 . 
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Figure 3: Sum of control coefficients as a function of total enzyme concentrations. Enzyme 
concentrations \Ej\ in /jM. For both enzymes SEi = 0.2[iM. 
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Figure 4: Sum of control coefficients as a function of total enzyme concentrations. Enzyme 
concentrations \Ei\ in /jM. For both enzymes SEi = —0.2/j,M. 
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Figure 5: Enzyme saturation as a function of total enzyme concentrations. Single enzyme 
saturation is measured as J / Eik cat ^y Double enzyme saturation is measured as {J/E\ks + 
J/E 2 k7)/2. Enzyme concentrations [Ei] in /j,M. 
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Figure 6: Measure of epistatic interaction effects on flux. Enzyme concentrations [E t ] in fiM. 
For both enzymes SEi = 0.02/iM. 
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Figure 7: Changes in control of a pathway due to alteration of catalytic turnover rate k cat 
(k C at(i) = 100 = 42, 000mm -1 ). Enzyme concentrations [Ei] in \xM. For both enzymes 8Ei = 
0.02^M. 
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Figure 8: Steady state substrate concentrations as a function of total enzyme concentrations 
in pathway pQ. Horizontal axes: (iM units. Vertical axes: mM units. ( see Appendix C ). 
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Figure 9: Steady state enzyme-substrate and free enzyme concentrations as a function of 
total enzyme concentrations in pathway pQ. Horizontal axes: fj,M units. Vertical axes: mM 
units ( see Appendix C ). 



40 




Figure 10: Failure of a modified summation conjecture.Enzyme concentrations [Ei] in fj,M 
( see Appendix C ). 



41 



REFERENCES 



Cornish-Bowden, A. (1987). Dominance is not inevitable. /. theor. Biol. 125, 333-338. 

Fontana, W. & Buss, L. W. (1996). The barrier of objects: from dynamical systems to 
bounded organizations. In Boundaries and Barriers, (Casti, J. & Karlquist, A., eds), 
pp. 56-116. Addison- Wesley, Reading, MA. 

Giersch, C. (1988). Control analysis of metabolic networks. Eur. J. Biochem. 174, 509- 
513. 

Grossniklaus, U., Madhusudhan, M. & Nanjundiah, V. (1996). Nonlinear enzyme kinetics 
can lead to high metabolic flux control coefficients: Implications for the evolution 
of dominance. /. theor. Biol. 182, 299-302. 

Heinrich, R. & Rapoport, T. A. (1974). A linear steady-state treatment of enzymatic 
chains. General properties, control and effector strength. Eur. J. Biochm. 42, 89-95. 

Kacser, H. (1987). Dominance not inevitable but very likely. J. theor. Biol. 126, 505-506. 

Kacser, H. (1991). A superior theory? J. theor. Biol. 149, 141-144. 

Kacser, H. (1995). The control of flux: 21 years on. Recent developments in Metabolic 
Control Analysis. Bioch. Soc. Trans. 23, 387-391. 

Kacser, H. & Burns, J. A. (1973). The control of flux. Symp. Soc. Exp. Biol. 27, 65-104. 

Kacser, H. & Burns, J. A. (1979). Molecular democracy: who shares the controls? 
Bichem. Soc. Trans. 7, 1149-1160. 

Kacser, H. & Burns, J. A. (1981). The molecular basis of dominance. Genetics, 97, 
639-666. 

Kacser, H., Burns, J. A. & Fell, D. (1995). The control of flux: 21 years on. The control 
of flux. Bioch. Soc. Trans. 23, 341-366. 

Keightley, P. D. (1996). A metabolic basis for dominance and recessivity. Genetics, 143, 
621-625. 



42 



Kholodenko, B., Rohwer, J., Cascante, M. & Westerhoff, H. (1998). Subtleties in control 
by metabolic channelling and enzyme organization. Mol. Cell. Biochem. 184, 311— 
320. 

Mayo, O. & Burger, R. (1997). Evolution of dominance: a theory whose time has passed? 
Biol. Rev. 72, 97-110. 

Porteous, J. W. (1996). Dominance-one hundred and fifteen years after mendel's paper. 
/. theor. Biol. 182, 223-232. 

Savageau, M. A. (1976). Biochemical Systems Analysis: A study of function and design 
in molecular biology. Addison-Wesley, Reading, MA. 

Savageau, M. A. (1992). Dominance according to metabolic control analysis: major 
achievement or house of cards? J. theor. Biol. 154, 131-136. 

Savageau, M. A. & Sorribas, A. (1989). Constraints among molecular and systemic 
properties: implications for physiological genetics. J. theor. Biol. 141, 93-115. 

Tralau, C, Greller, G., Pajatsch, M., Boos, W. & Bohl, E. (2000). Mathematical treatment 
of transport data of bacterial transport systems to estimate limitation in diffusion 
through the outer membrane. /. theor. Biol. 207, 1-14. 



43 



